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ABSTRACT 

As part of a radial velocity survey of low Galactic latitude structures that we under¬ 
took with the 2dF spectrograph on the AAT, we present the radial velocities of more 
than 1500 Red Giant Branch and Red Glump stars towards the centre of the Ganis 
Major dwarf galaxy. With a mean velocity of 72 ± 7kms“^ at a Heliocentric distance 
of 5.5 kpc and 114±2 kms“^ at 8.5 kpc, these stars present a peculiar distance - radial 
velocity relation that is unlike that expected from thin or thick disc stars. Moreover, 
they belong to a kinematically cold population with an intrinsic dispersion that may 
be as low as llj'j^kms”^. A comparison of the velocity distribution obtained in this 
work with previous studies, shows the importance of using our new reduction pipeline 
and averaging the velocities obtained from different templates. 

The radial velocity distribution is used to select Canis Major stars in the UCAC2.0 
proper motion catalogue and derive proper motions in Galactic coordinates of 
{ftijfib) = (—3.6 ± 0.8mas/yr, 1.5 ± 0.4mas/yr) for the dwarf galaxy, which af¬ 
ter correcting for the reflex solar motion along this line-of-sight gives (/i(,p,(,) = 
(—6.8±0.8mas/yr, 0.8±0.4 mas/yr), corresponding to a prograde orbit with a tangen¬ 
tial velocity of ~ 235kms“^ at the average distance of ~ 7.2 kpc. All these kinematic 
constraints can be reproduced in simulations of the accretion of a dwarf onto the 
Galactic disc. Such a process could also be responsible for the Monoceros Ring that 
has recently been shown to encompass the Galactic disc. However, without constraints 
on the kinematics of the tidal arms emerging from the Ganis Major dwarf, it is not 
yet possible to definitively prove a link between the two structures. 

Key words: Galaxy: structure - galaxies: interactions - galaxy: individual (Ganis 
Major) 


1 INTRODUCTION 


The advent of all sky surveys is revealing numerous struc¬ 
tures toward the edge of the Galactic disc. The SDSS re¬ 
vealed a ring-like structure, the so-called Monoceros ring 
(Mon ring), that was later shown to encompass part of the 
Galactic disc (e.g . |Newberg_^_aL| 200jjYamw 1200,^ 


llbata et al.ll200^ : lcranee'^1 .112003: IConn et aI Jl200,5^ . Sim¬ 

ilarly, the 2MASS infrared catalogue was used to unveil the 
existence of a diffuse structure in the d irection of the Tri- 

S ’um and Andromeda constellations iRocha-Pinto etld] 
and to reveal the presence of a dwarf galaxy below the 


120041) and to reveal tne presen ce ot a dwart galaxy below tne 
disc in the CMa constellation iMartin et al.l2004al . hereafter 
Paper I). 


This latter structure was first identified as an overden¬ 
sity of Red Giant Branch (RGB) stars located just below 
the Galactic disc at {l,b) (240°,—8°). The analysis of 

archival Colour-Magnitude Diagrams (CMD) of open clus¬ 
ters that are fortuitously located around this region con¬ 
firmed thi s stellar overdensity ex tends over 20° in Galactic 
longitude llBellazzini et aBl2004l . hereafter Paper II) while 
main sequence fitting revealed this population has an inter¬ 
mediate age (4-10 Gyr). Applying a Tip of the Red Giant 
Branc h algorithm yielded a distance of 7.2 ± 0.3 kpc for the 
dwarf llMartin et al.ll2f)04ll hereafter Paper III), making it 
the closest galaxy from the Sun. Interestingly, Canis Major 
lies at a comparable Galactocentric distance to the Sagit¬ 
tarius dwarf. 
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As an alternative interpretation, iMomanv et alJ (l2004l 
explain this stellar structure as a signature of the Galactic 
warp. Using the UCAC2.0 proper motion catalogue, they 
derive a proper motion that, within sizable uncertainties, is 
compatible with a Galactic population on a prograde disc¬ 
like orbit: (pa cos{S),ns) = (—1.7±2mas/yr, 3.1±2 mas/yr). 
However, comparison of deep optical photometry CMDs 
with Galactic models reveals the structure is incompati¬ 
ble with known Galactic structure (Paper II) and yields a 
narrow e xtent on the line of sight, wit h a FWHM of only 
1.92 kpc JMartmez-Delgado et al.ll200R^ . difficult to achieve 
with a conventional warp. Moreover, the radial velocities of 
a small sample of RGB stars in CMa show they belong to 
a kinematically cold population that, if on a circular orbit, 
would rotate around the Milky Way at a rotational velocity 
of ~ 160kms“^ that is low for disc stars (Paper III). 

The existence of a dwarf so close to the Galactic disc, 
along with structures like the Mon Ring and the TriAnd 
feature raises questions on their role in the formation of 
the Galactic (thick) disc (Paper I). However, answering this 
question hrst requires us to determine if there i s a lin k be¬ 
tween these three structures. IPefiarrubia et all ^2^^!^^ used 
all known distance, radial velocity and proper motion data 
to constrain simulations of the formation of the Mon Ring 
by an accretion event. Interestingly, they predict the pro¬ 
genitor of the Ring should be located around the position of 
CMa {I ~ 245° and b ~ —18°) but at twice the estimated 
distance of the CMa dwarf. Therefore, it remains unclear 
whether the two structures are linked or not. 

This paper is the hrst in a series that aims at analyzing 
the kinematics of these different low latitude Galactic struc¬ 
tures. In particular, we aim at constraining the orbit of the 
accreted dwarf and determining if the Mon Ring could be a 
by-product of this accretion. Here, we present our complete 
sample of radial velocities of Red Giant Branch and Red 
Clump stars in the CMa overdensity and show it presents 
non-Galactic features. In section 2 we discuss the data set 
we use for our radial velocity study of section 3. These radial 
velocities are used in section 4, to select CMa stars from the 
UCAC2.0 catalogue and determine the proper motion of the 
dwarf. The obtained kinematics are then used in section 5 
to constrain simulations of the accretion of a dwarf onto the 
Galactic disc. Conclusions are presented in section 6. 

In the following, all the J, H, K magnitudes from 
2MAS S have been corrected from ex ti nction using the maps 
from ISchlegel. Finkbeiner fc D ayis| jl998|). modihed by 
the as ymptotic correction from iBonifacio. Monai fc BeersI 
i200(tl . We also assume that the Solar radius is Rq = 8 kpc, 
that the LSR circular velocity is 220kms“^, and that the 
peculiar motion of the Sun is ( Up = 10.00 kms~^, Vo = 
5.25kms“^,Wo = 7.17kms~^: IPehnen fc BinnevlflQ^ . 
Except when stated otherwise, the radial velocities are the 
observed Heliocentric radial velocity, not corrected for the 
motion of the Sun. 


2 OBSERVATIONS AND REDUCTION 

The aim of the low Galactic latitude survey we undertook 
with the 2dF spectrograph was to determine the kinematics 
of the different asymmetries that appear in the distribution 
of Red Giant Branch (RGB) stars (Figures 3 to 5 of Pa- 


Canis Major fields reference field 




Figure 1. Map of the target Red Giant Branch stars for the 
CMa region (left panel) and the reference field (right panel). All 
RGB stars of these regions are shown as dots and target stars 
are circled. Extinction contours for E{B — V) = 0.1, 0.25 and 
0.55 are also shown as thick lines. Target Red Clump stars follow 
the distribution of target RGB stars but are not shown to avoid 
overcrowding the plots. 

per I) and in particular kinematics of the Ganis Major dwarf 
galaxy. Hence, four fields out of the ~ 15 in our survey were 
directly aimed at the centre of the CMa overdensity as found 
in Paper I: {l,h) = (240.0°,-4.8°), {l,h) = (240.0°,-6.8°), 
{l,b) = (240.0°,-8.8°) and (Z,6) = (240.0°,-10.8°). We also 
observed another field on the opposite side of the Galactic 
disc, at {l,b) = (240.0°,-1-8.8°) to provide a reference to 
compare with the CMa fields. 

In each two-degree field, our primary targets were the 
RGB stars selected from sample A of Paper I and within 4 
to 20 kiloparsecs from the Sun^. These stars are highlighted 
on the map of sample A stars shown on Figure 1 and mainly 
fall in regions of reasonable extinction {E{B — V) < 0.4). 
Only the field closest to the Galactic disc with —5.8° < 
b < —3.8° suffers from significant variable extinction. The 
target stars are also shown on the CMD of each field in 
Figure 2. As the CMa structure has a prominent Red Clump 
(RC, see Papers I and H), we placed the remaining fibres to 
observe a sample of RC stars at the estimated distance of 
CMa: 5 kpc < Dq < 8 kpc, assuming a distance of 7.2 kpc 
for the bulk of the overdensity (Paper HI) and an absolute 
magnitude of Mk = —1.5 ± 0.2 for the RC population of 
CMa (Paper H), independent of colour which was chosen 
within the range 0.5 < J — K < 0.65. The corresponding 
selection box has also been drawn on Figure 2. It can be 
seen that even in fields with increasing extinction values 
(panels c and d), the selection box remains centred on the 
Red Clump zone, hinting at a good extinction correction for 
stars at this distance. 

The observations were obtained during the nights of 7- 
12 April 2004. We employed two different spectrograph set- 

^ We recall that the distance to RGB stars weredetermined_^ 
th^photometric parallax technique presented by |Majewsld^ alJ 
l|200^ that we applied to a sample of CMa stars in Paper III. 
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Table 1. Templates used for cross-correlation of first spectro¬ 
graph observations. 


Figure 2. 2MASS infrared Colour Magnitude Diagrams of the 
four target fields around Canis Major {a to d panels) and the 
symmetric field above the Galactic disc (panel e). On each panel, 
the observed RGB stars are shown as circled dots and the selection 
box from which Red Clump stars were chosen is drawn in white. 


tings, with the 1200V grating on spectrograph 1 (covering 
4600-5600 A at 1 A/pixel) and with the 1200R grating on 

o o 

spectrograph 2 (covering 8000-9000 A, also at 1 A/pixel). 
The observations have been reduced usin g the 2dF Data Re - 
duction package provided by the AAO IITavlor et alJll99^ 
for the correction from the flat fields, the extraction of the 
spectrum of each fibre and the sky subtraction. 


2.1 Reduction of first spectrograph spectra 

We use a custom-made reduction pipeline to correct the 
asymmetry of the Line Spread Function (LSF) of the 2dF 
with the first spect rograph setti n gs. T his pipeline is de¬ 
scribed in detail in iMartin et al.l l^200!^^ and can be sum¬ 
marized as follows: 

(i) the asymmetry of the LSF is modeled across the 
CCD for each observed field; 

(ii) the spectrum of a given fibre is calibrated using the 
corresponding LSF model; 

(hi) for each observed spectrum, the model is also used 
to generate template spectra. Hence, when doing a Fourier 
cross-correlation to determine the radial velocity of an ob¬ 
served star, both the observed spectrum and the templates 
are deformed in an identical way. 

We have shown that this procedure ensures that we avoid 
systematic offsets higher than ±5kms“^ due to deforma¬ 
tions of the LSF. 

Each observed spectrum was cross-correlated with up 
to 9 different templates of giant stars, from K4III to 
M4III spectral types, generated from high res olution spectra 
from the UVES Paranal Obs ervatory Project iBagnulo et alJ 
l200,8l . see lMartin et aliboOi'il and Table 1 for the list of tem¬ 
plates). RGB stars are cross-correlated with all nine tem¬ 
plates and have typical nncertainties of ~ 4 — 7kms“^ 


Star 

spectral type 

radial velocity (kms ^)“ 

HD 145206 

K4HI 

-46.0 

HD 167818 

K5HI 

-16.9 

HD 149447 

K6HI 

-2.1 

HD 89736 

K7HI 

16.0 

HD 92305 

MOIII 

-22.4 

HD 102212 

Mini 

50.7 

HD 120052 

M2HI 

64.2 

HD 224935 

M3HI 

-11.8 

HD 11695 

M4HI 

1.5 


“ Obtained using the SIMBAD database, operated at CDS, Stras¬ 
bourg, France. 

The high resolution spectra of these stars are extracted fr om the 
UVES Paranal Observatory Project iBaenulo et al 1I2 oo3i 

for each derived radial velocity while RC stars are cross- 
correlated with the first five templates (from K4III to MOIII 
spectral type), since their spectra are unlikely to correspond 
to a higher spectral type, and have typical uncertainties of 
~ 8 — 15kms“^ for each template. Instead of choosing one 
of the template-specific derived radial velocities, we use a 
weighted average of the nine values to avoid any systematic 
offset that could be due to particular features in one of the 
templates. Hence, the Heliocentric radial velocity, Vr, of an 
observed star is given by: 


nb temp 


nb temp 




( 1 ) 


where Vr,i is the radial velocity derived from the cross¬ 
correlation of the observed spectrum with the ith template 
and (Ji is the related uncertainty, given by the iraf function 
f xcor. To judge the homogeneity of the different Vr,i around 
Vr, we also calculate the weighted dispersion, ct(,, of the Vr^i 
around Vr'- 


nb temp 




nb temp - 

--- {Vr-Vr,i) 


( 2 ) 


Stars with poorly determined radial velocities were 
eliminated by keeping only stars with ct' < 5kms“^ which 
represent more than 90 percents of our sample. 


2.2 Reduction of second spectrograph spectra 

Only RC stars were observe d with the second sp ectrograph 
settings. It has been shown in iMartin et alJ (l200?t) that these 
settings produce no important deformation of the LSF and 
do not require the application of the reduction pipeline we 
used for spectrograph 1 data. For each star, we fit a Gaussian 
model for each of the three lines of the Ca II triplet and 
derive a velocity. We use the weighted average of these three 
velocity values as the radial velocity of the star and compute 
the weighted dispersion ct),. Stars with ct), < 6.0kms“^ are 
kept as valid stars for the second spectrograph settings (once 
again this represents more than 90 percent of the sample). 

On both spectrographs, and to account for systematic 
offsets of no higher than (T 2 dF = 5kms“^ that may remain 
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Figure 3. The top panel shows the radial velocity distribution of 
RGB stars near the centre of CMa (5 kpc < Dq < 11 kpc). A peak 
of stars is present at ~ 115kms“^. The lower panel shows the 
position of these RGB stars in phase space, with the typical error 
on both radial velocity and distance reported on the left. The 
peak is produced by a group of stars clustered between 7 kpc and 
9 kpc. The expected position of a population orbiting the Milky 
Way in a circular orbit at Urot = 240kms~^ (full line), Urot = 
200kms~^ (dashed line) and Urot = 160kms~^ (dashed-dotted 
line) has been overplotted for comparison. The stars selected to 
produce the histogram of the upper panel are those between the 
two dotted line. 

after applying our reduction pipeline, we increase the uncer¬ 
tainty cr(, of each star to reach the total uncertainty <7v on 
each radial velocities: 

(Jv = \/ o'^ + cr|jp (3) 

Finally, we note that the observed spectra are of suffi¬ 
ciently good quality to avoid significant variation of cr„ with 
the magnitude of the target stars ((j„ varies between 6 and 
8kms“^ over the 7 < < 13 range of our sample). 

3 KINEMATIC RESULTS 

3.1 The properties of the RGB sample at the 
centre of CMa 

Our sample of RGB stars show no significant change with 
Galactic latitude hence, we merge the data of all the fields 
with b < 0°. The distribution of radial velocities of the 228 
RGB stars in these four GMa fields is shown on the top panel 


Paper III templates This paper's templates 



Figure 4. Distribution of the radial velocities produced by the 
different template of Paper III (left panel) and this paper (right 
panel) for the same RGB star. The thick lines represent the radial 
velocity used in the two papers and the dotted line represent the 
ilcTi, limit. The reduction pipeline used here yields much more 
clustered values than before. See the text for more details. 


of Figure 3 selected around the estimated distance to the 
structure (5 kpc < Dq < 11 kpc). The distribution is rather 
broad, but has a well-defined peak of stars at ^ llSkrns”^. 

3.1.1 Comparison with the \Martin et al\ l200Jll) results 

In Paper III, we presented a similar analysis of the ra¬ 
dial velocity of RGB stars near the centre of Ganis Ma¬ 
jor based on only the two central fields of the present sur¬ 
vey [(Z,6) = (240.0°,-6.8°) and (Z, 6) = (240.0°,-8.8°)]. 
However, in that analysis the distribution of radial velocity 
showed two narrow peaks, centred on 60 and llOkms”^, 
that are not reproduced here. The second peak is visible 
on Figure 3, but the first one has completely disappeared. 
This peak was in fact artificially produced by template is¬ 
sues resulting from a fluctuating LSF asymmetry during the 
different observation nights. 

This is illustrated on Figure 4 where, for one of the stars 
belonging to the first peak of Paper III, we present the radial 
velocities derived from the different observed templates (left 
panel) and those derived using the artificial templates to 
correct for the LSF asymmetry (right panel). Using our new 
reduction pipeline, most of the derived radial velocities are 
well clustered around Vr as defined in equation (1). Only two 
templates are over the 1-a limit but they have high radial 
velocity uncertainties and hence, are not significant in the 
determination of the combined Vr. On the other hand, the 
radial velocities derived for the different templates used in 
Paper III are widely scattered (over 35kms“'^), and are not 
uniformly distributed. In particular, the two low uncertainty 
values are clustered around ^ 60kms“'^ while all the other 
values are clustered not far from the new radial velocity 
derived by the pipeline. Since in Paper III we only used the 
radial velocity that had the lowest uncertainty as the radial 
velocity of the studied star, this star (as well as multiple 
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Figure 5. The radial velocity distribution of those stars within 
7.5 and 11 kpc from the Sun in the four CMa fields. The peak in 
the distribution is well reproduced by a Gaussian model centred 
on 114±2 km s“^ and with an intrinsic dispersion of ll^j kms“^. 


other ones) ended up artificially populating a peak of stars 
at ~ 60kms“^. 


3.1.2 The Canis Major radial velocity signature 

Looking at the position of the RGB stars in phase space 
(bottom panel of Figure 3) reveals the peak in the distribu¬ 
tion is produced by a clump of stars with 7.5 kpc < Dq < 
11.0 kpc. A more detailed analysis of this distance range 
(Figure 5) reveals a population with a narrow dispersion. 
Comparing the data between 85kms“^ and 145kms“^ with 
a Gaussian model using a maximum likelihood technique re¬ 
veals this population is centred on 114 ± 2kms“^ and has 
an intrinsic dispersion of lljijkms”'^ when accounting for 
radial velocity uncertainties as in equation (3). As we ex¬ 
plained previously in Paper III, the radial velocity of this 
population is not easily compatible with a disc-like popula¬ 
tion since it would be orbiting the Galaxy with a rotational 
velocity of only ~ 160kms“^. Given the distance and pe¬ 
culiar velocity of this population, we assume that it corre¬ 
sponds to the CMa dwarf. Changing the lower distance cut 
in the range 7 kpc to 8 kpc does not substantially modify 
these results (position and dispersion change by less than 
2 kms“^). 

At closer distance, even though 2MASS starcounts of 
RGB stars betray the presence of the CMa dwarf (Paper III), 
most of the stars follow a more disc-like velocity. Yet, a group 
of these stars seems aligned in phase space, from Dq 6 kpc 
and Vr ~ 70kms“^ to the position of the clump of stars at 
higher distance. The low number of RGB stars in our sample 
prevents any firm conclusion on the reality of this feature, 
even though it seems to have a low dispersion. In the same 
distance range, a more diffuse group of stars also appears at 
higher radial velocity {vr > 140 kms“^). Given the relatively 
sparse nature of these stars they could simply be the tail of 
the distribution of disc stars and/or could possibly represent 
a contamination from halo stars along these sight lines. 

Another interesting feature of the phase space diagram 
is the group of stars at high distances (mainly 11 kpc < 
Dq < 15 kpc corresponding to 16.5 kpc < Doc < 20 kpc) 
that seems to be disconnected from the CMa population 


Vr (km/s) (km/s) v^. (km/s) 

0 50 100 150 0 50 100 150 0 50 100150 




Vr (km/s) 

Figure 6. Phase space distribution of stars in the RC sample 
(bottom panel). A significant part of the stars follows a distance 
- radial velocity relation that is confirmed when analyzing the 
radial velocity distribution of stars in 1 kpc bins (histograms of 
top panels, with the fits corresponding to the values of Table 2). 
The positions of this peculiar population have been overplotted 
as filled circles in the bottom panel, with error bars representing 
1(7 uncertainties on these values. 

by a gap at Dq ~ 11 kpc. The Galact ic disc is known to 
possess a cut off at Dec ~ 15 kpc (e.g. iR.unhv et alJlTflO^ 
and hence, no disc stars are expected at the distance of this 
group. Moreover, given the high radial velocity of these stars, 
it seems unlikely that this group is composed of misidenti- 
fied dwarfs. In fact, since these stars appear to be located 
between 16.5 to 20 kpc from the Galactic centre we believe 
we have uncovered part of the Monoceros Ring behind the 
CMa dwarf. The radial velocity signature of this population 
is compatible with previous work but we defer a more thor¬ 
ough analysis of the presence of the Monoceros Ring in all 
our low latitude fields to another contribution (Conn et al. 
in preparation). 

3.2 The RC sample at the centre of CMa 

Most of the RC stars were observed in the {l,b) = 

(240.0°,—8.8°) field and should give a more precise view 
of the kinematics at the core of the CMa dwarf. While re¬ 
stricted to a shorter range in distance, the RC sample has 
the advantage of containing ~ 1350 stars at the distance of 
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Table 2. Parameters of the Red Clump radial velocity distribu¬ 
tion fits. 


Distance bin 
(kpc) 

< 'Cr,CMa > 

(kms“^) 

f^CMa 

(kms“^) 

’^r,cont ^ 

(km s~^) 

^cont 

(km s~^) 

5 < Do < 6 

72 ±7 

24 ±4 

75 ±6 

42 ±4 

6 < Do < 7 

87 ±4 

20 ±3 

76 ±5 

54 ±4 

7 < Do < 8 

102 ±3 

16 ± 3 

71 ±5 

46 ±3 


Column 1 and 2 state, respectively, the mean velocity and the 
intrinsic dispersion of the CMa population. Columns 3 and 4 list 
the same parameters for the contaminating population. 

the group of RGB stars that seem to be aligned in phase 
space. The distribution of these RC stars in phase space 
(bottom panel of Figure 6) presents a behaviour that is very 
similar to the RGB sample, with an overdensity of stars 
that follow a distance - radial velocity relation. The distri¬ 
butions of the radial velocity of these stars for 1 kpc bins are 
shown on the top panels of Figure 6. We suspect two popu¬ 
lations are present in each bin: contaminating stars that fall 
in our GMD selection box of RG stars and genuine Canis 
Major RC stars. To determine whether it is worth fitting 
a double Gaussian mo del to the distribution , we u se the 
KMM test presented in lAshma.n. Bird fc Zenll The 

probabilities P that the data are better represented by a 
single Gaussian instead of two are P(5 < Dq < 6) = 1.0, 
P{6 < Dq < 7) = 6 ■ 10“® and P(7 < Dq < 8) < 10“^ for 
the three distance bins. The last two bins are hence inconsis¬ 
tent with a single Gaussian model and only the closest one 
would be better characterized by a single Gaussian. How¬ 
ever, since the distance - radial velocity relation still seems 
to be present in this bin, we believe that the two populations 
appear in this bin but that their mean velocity is located at 
around the same position. Hence, we proceed in fitting a 
double Gaussian model to the data in all bins, knowing that 
uncertainties on the derived values for the first bin will be 
large. 

We use a maximum likelihood technique to fit the dou¬ 
ble Gaussian model. The uncertainties (j„ on each radial ve¬ 
locity were taken into account whereas no uncertainty was 
assumed for RC distance since the only important source of 
error comes from the absolute magnitude of the population 
and a change in this value would only shift the sample as a 
whole. Since the two populations are most clearly separated 
in the 7 < Dq < 8 kpc, we use the stars in this bin to de¬ 
termine the proportion of RC stars that belong to the Canis 
Major population. The fit yields to a 40 percent proportion 
of CMa stars, well within the 25 to 50 percent determined 
bv lBella.zzini et !iD (Emi). Since we expect this proportion 
to be roughly constant over our 3 kpc sample (see Figure 2 of 
Paper H), we also adopt this value for the two closest bins. 
The four parameters of the fits - mean velocity < Ur.CMa > 
and intrinsic dispersion acMa of the Canis Major population 
and mean velocity < Ur,cont > and intrinsic dispersion CTcont 
of the contaminating population - are presented in Table 2. 

Each distance bin contains a low dispersion population 
with an increasing mean radial velocity and a broad popula¬ 
tion of constant mean velocity. This latter population shows 
a velocity dispersion (Tcont comparable to Solar neighbour¬ 


hood velocity dispersions of the thick disc (au,av,aw) = 
(63±6, 39±4, 39±4) kms”*^ dSoubiran. Bienavme fc SiebertI 
|2003l- We believe these stars are disc stars and dwarfs close 
to the Sun that fall in the RC region of the (J-K,K) GMD 
and contaminate our sample. On the other hand, the other 
population that we believe to be composed of Canis Major 
stars is kinematically cold, with an intrinsic internal disper¬ 
sion around 16 ± 3kms“^. It is likely that the increase in 
ccMa in the other distance bins is only due to difficulties for 
the fit to disentangle the two populations. The other notable 
feature of these stars is their increasing radial velocity with 
distance, from 72±7 kms“^ to 102T3 kms”*^ over only 3 kpc. 
Comparison with the group of RGB stars we identified as 
Canis Major stars at higher distance reveals striking simi¬ 
larities with both RC and RGB samples being kinematically 
cold and having compatible intrinsic dispersion. Moreover, 
the distance - radial velocity trend that appears in the RC 
sample is nicely prolonged to higher distances by the group 
of Canis Major RGB stars at Dq ~ 8.5 kpc and overlap 
with the group of RGB stars that seem to be align in phase 
space for the Dq > 6 kpc region. It is unlikely that a disc¬ 
like population could show characteristics that are similar to 
these stars given the high shift in radial velocity over only 
a few kiloparsecs. Indeed, if the considered population was 
rotating around the Milky Way on a circular orbit, it would 
have a rotational velocity as high as Vrot ^ 220kms“^ at 
5.5kpc and as low as Vrot ^ 160kms“^ at 8.5kpc. On the 
other hand, correlation in phase space and low dispersion 
in radial velocity, two characteristics that the CMa popu- 
lation shows, are typical of an accretion process (see e.g. 
Ilbata et al ] ll997l. for the Sgr dwarf). 

To further test this conclusion, we compare the o bser- 
vations with the Besangon model jR.obin et al 1 1200311 . We 
simulate the J,H,K Colour-Magnitude Diagram of the re¬ 
gion 239° < I < 241° and -4° < b < -12°, from which we 
extract Red Clump stars with the same colour cuts as those 
we used to define our RC sampie. We then calculate the dis¬ 
tance to the stars in the model as we do for the observed 
stars, assuming an absolute magnitude of Mk = —1.5. Of 
course, since these stars do not belong to the CMa dwarf, 
their estimated distance is different from their distance pro¬ 
vided by the model. From the ~ 150 stars in the model for 
the selected region, four fifths are in fact located at less than 
1 kpc from the Sun. The distribution of these stars in phase 
space (Figure 7) therefore represents that of the close disc 
stars that contaminate our sample of CMa RC stars. They 
do not follow the distance - radial velocity relation that we 
attributed to the CMa dwarf and are, on the contrary, lo¬ 
cated at a lower velocity. A similar behaviour is observed for 
RGB stars in the model. 

3.3 The {l,b) = (240.0°,-h8.8°) field 

Another check of the peculiar radial velocity of the stars in 
the CMa structure is given by the field that is symmetric to 
the CMa fields on the other side of the Galactic disc and that 
does not show the features that appear at the centre of the 
CMa dwarf (Figure 8). The RGB sample for this field con¬ 
tains stars that are mostly located at Dq < 6 kpc with the 
expected radial velocity of a population orbiting the Milky 
Way at Vrot ~ 200kms“^. Only 11 stars fall within the dis¬ 
tance cut that was used to produce Figure 5 near the centre 
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Figure 7. Phase space distribution of stars in the Besangon 
model that fall in our selection box for RC stars in the region 
of the observed 2dF fields. For the top panel, the distance of the 
stars in the model was used and shows that most of the contam¬ 
ination in the observed sample is due to stars close to the Sun. 
For the bottom panel, we calculate the distance as we do for the 
observed sample of RC stars. A direct comparison of this panel 
with Figure 6 shows that the distance - radial velocity relation 
in our sample is not predicted by the model. 



■^50 0 50 100 150 200 


Vr (km/s) 

Figure 8. Phase space distribution of RC (top panel) and RGB 
stars (bottom panel) in the {l,b) = (240.0°,+8.8°) field. Neither 
of the two show the features that are present in fields centred on 
the CMa dwarf galaxy. 


of CMa. Moreover, the remaining stars at higher distances 
show no sign of clustering. The RC stars in this field fol¬ 
low a broad distribution with a dispersion of ~ 50kms“^, 
centred around ~ 60kms“^ over the 3kpc sample, similar 
to what we identify as the contaminating population in the 
CMa field, and also similar to what would be expected from 
close stars artificially spread over the 3kpc distance range. 
Therefore, the peculiar features that appear at the centre of 
CMa have to be due to the CMa population. 


4 UCAC2.0 PROPER MOTIONS OF THE CMA 

DWARF 

iMoman^^^^J j200^ used the UCAC2.0 catalogue 
IZaciiaria^^ to show that the CMa overdensity 

of RGB stars rotates around the Galaxy in a prograde mo¬ 
tion. However, given the sizable uncertainties on the proper 
motion values of individual stars and the contamination of 
disc stars in the sample, they were only able to give a rough 
estimate of the proper motions that convert in Galactic co¬ 
ordinates as: m = —3.5mas/yr and fib = —O.lmas/yr with 
uncertainties of ~ 2.0mas/yr. 

With the addition of radial velocities of stars in the di¬ 
rection of CMa, we are now in a position to define a less 
contaminated sample of stars that belong to the dwarf. We 


first determine the best linear fit to the distance - radial 
velocity relation of the RC sample by using the three radial 
velocity measurements for the three distance bins. We then 
extract from UCAC2.0 the proper motions of the RC stars 
that are located within TlOkrns”^ of this linear fit. This 
selection should ensure that a sizable proportion of our ob¬ 
jects indeed belong to the CMa dwarf. To account for the 
distance Dq of the different stars, we normalize each proper 
motion fivcAC as if the star was located at the Heliocentric 
distance of the CMa dwarf (DQ^cMa = 7.2 kpc): 


n ~ nvcAC ■ 


Dq 

DQ,CMa 


(4) 


The distribution of these fi proper motions is shown on the 
bottom left panel of Figure 9 and is consistent with the 
iMomanv et alJ ll2nn4h measurements. 

We use a maximum likelihood technique to deter¬ 
mine the best fit of this distribution by a two di¬ 
mensional Gaussian function. The best mean value is 
(jlijjlb) = (—2.5 mas/yr, 1.5 mas/yr) with internal disper¬ 
sions of (cr^ijcr^b) = (4.8mas/yr,3.7mas/yr). To analyze in 
more depth the proper motions of the CMa population, we 
use these values to eliminate outliers and study the fii dis¬ 
tribution of stars within cr^f, of Jib (top panel of Figure 9) 
and the fib distribution of stars within of Jli (right panel 
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Figure 9. The bottom left panel shows the distribution in proper 
motion space of CMa RC stars selected along the distance - radial 
velocity feature of Figure 6. A black pixel corresponds to 8 stars. 
The centre of the distribution (assuming Gaussian distributions) 
is overplotted as a dot. The distribution of those stars within 1 a 
of this point in the p(, (respectively p() direction where used to 
produce the distribution of CMa stars in p; (resp. on the top 
(resp. right) panel. The mean position derived from the Gaussian 
fits of the whole distribution are shown as arrows. A wing in the 
distribution in p; betrays the presence of two populations that 
we fit by a double Gaussian model. 

of Figure 9). This latter distribution is well centred around 
Jib (the arrow) and the value 1.5 ± 0.5mas/yr can be taken 
as the proper motion of the CMa dwarf in Galactic lati¬ 
tude. On the other hand, the distribution of proper motions 
in Galactic longitude is not centred around Jli. A peak ap¬ 
pears at lower proper motion, with a wing that extends to 
positive values. Considering that, even with our cut in ra¬ 
dial velocity, contaminating stars from the disc should still 
be in our sample of stars, we use a maximum likelihood 
technique to fit a double Gaussian model to the distribu¬ 
tion (the KMM test yields a low probability of 10”^ that 
the population is better represented by a single Gaussian). 
This reveals two distinct populations are present: similar to 
what was observed in radial velocity, one has a large dis¬ 
persion (4.7mas/yr), is centred on /i; = —1.9mas/yr and 
accounts for one fourth of the total number of stars while 
the other has a much narrower dispersion (1.5 ± 0.4 mas/yr) 
and is centred on fj.i = —3.6 ± 0.8 mas/yr. Since one expects 
the contaminating stars in the sample to be at different dis¬ 
tances and hence have different proper motions, we believe 
the broad population corresponds to these stars; so the nar¬ 
row population has to be due to the CMa dwarf. 

For comparison, we show on Figure 10 the distribution 
of RC stars in our sample that are more than ±20kms“^ 
away from the linear fit of the Canis Major radial velocity 
and that should correspond to the contaminating popula¬ 
tion in the sample. This population has a proper motion 
of (tti,cont,Ms>,cont) ~ (—2.0 mas/yr, 2.0 mas/yr) offset from 



(mas/yr) 

Figure 10. The contours show the distribution of RC stars in 
our sample that are more the 20kms~^ away from the distance 
— radial velocity relation followed by CMa stars. This population 
is centred around l^b,cont) ~ (-2.0mas/yr, 2.0mas/yr) 

and offset from the CMa population whose mean proper motion 
is represented by the black dot. 


the CMa proper motion, especially for the proper motion in 
Galactic longitude. 

The derived proper motion values for the CMa dwarf 
{ni,nb) = (—3.6 T 0.8mas/yr, 1. 5 ± 0.4 mas/yr) are com¬ 
patible with the iMomanv et al.l (|2003) values but with 
much lower uncertainties. Correcting from the reflex solar 
motion along the line of sight yields {fi'i, Fb) — (~6.8 ± 
0.8 mas/yr, 0.8 ± 0.4 mas/yr), corresponding to a tangential 
velocity of 234kms“^ at the mean distance of 7.2 kpc. 


5 SIMULATIONS 

The kinematic information on the CMa dwarf provides in¬ 
teresting constraints on the accretion process that the dwarf 
is undergoing. In this section, we revisit the simulations of 
Paper I, this time using only constraints that can be directly 
linked to CMa. Indeed, since there is no definite proof of a 
link between CMa and the Mon Ring, we choose not to use 
the numerous kinematic and/or positional data on the Mon 
Ring for our simulations^. 

The simulations are performed in the same way as 


^ Moreover, the existence of the Mon Ring behind CMa seems to 
hint against a direct link between the two structures and it would 
be the tidal arms of the dwarf, wrapped a few times around the 
Galaxy, that would prod uce the ring-like structu re. If such a sce¬ 
nario is plausible (see e.g.|Penarrubia^_et_^lj200^, the constraints 
provided by the Ring data on the accretion parameters of the 
dwarf are weak. 
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Figure 11. Our best simulation of the accretion of the CMa dwarf 
viewed from the North Galactic Pole. The Sun is represented by 
a dotted circle and the three dashed circles represent distances 
of 10, 15 and 20kpc from the Galactic centre. The main body of 
the GMa dwarf is visible at (X, Y) ~ (—12 kpc, —6 kpc) and the 
tidal arms produced by the accretion are wrapped multiple times 
around the Milky Way. 

in iMartin et al.l (|20o3); using the fast an d momentum- 
conse rving tree code integrator falcON dPehnerJ 1200(1 
|20 o 3)- The simulation that best fits the kinematic data at 
the centre of CMa is shown from the North Galactic Pole on 
Figure 11. It is produced by a dwarf galaxy modeled by a 
King model of 5 x 10® Mq with a tidal radius of 2.7 kpc and 
Wo = 3.25 and that is accreted onto the Milky Way dur¬ 
ing ~ 3 Gyr. The simulation reproduces at the same time 
the overdensity of stars that revealed the dwarf in 2MASS, 
the proper motions we measured in this paper (with val¬ 
ues of Pi ~ —4mas/yr and p(, ~ 2mas/yr) and, above all, 
the distance - radial velocity gradient that appears in RC 
(and possibly RGB) stars. This is shown on Figure 12, with 
the values measured in this paper plotted as squares and 
triangle on top of the particles of the simulation in the re¬ 
gion of our 2dF fields. Moreover, the velocity dispersion of 
the particles reproduces that measured for the RGB sample 
(llkrns”^). Of course, the parameters of the initial dwarf 
are certainly not the only ones that can reproduce the ra¬ 
dial velocities at the centre of the dwarf but it should be 
noticed that the simulated position and velocities are very 
near what is observed. 

The group of particles that is visible at higher dis¬ 
tance (9 kpc < D < 12 kpc) is produced by one of the 
tidal arms of the dwarf that is wrapped a few times around 
the Milky Way. Even if the group is not exactly at the 
same distance as what we identify as the Mon ring behind 
the GMa dwarf, it is tempting to explain the two obser¬ 
vations by the same accretion process. It is also interest¬ 
ing to notice that the average location of these tidal arms 
in Figure 11 are close to the distance range of the Mon 



Figure 12. Distribution in space phase of the particles (dots) of 
our best simulation at the location of the 2dF fields we present in 
this work (239° < I < 241° and —12° < b < —4°). The observed 
values of radial velocity are shown as squares for values derived 
from the RC sample and a triangle for the one derived from the 
RGB sample. They are well matched by the particles. The group 
of star at higher distances is produced by the multiple wraps of 
the tidal arms of the dwarf around the Milky Way. 

ring (10 kpc < Dec < 18 kpc in the anticentre direction) 
while the radial velocity of the arm reveals it is only shifted 
by ~ 30kms~^ co mpared to the SDSS radial velocities of 
lYa.nnv et a,1.l The more distant portions of the arms, 

that extend at 20-30 kpc from the Galactic centre could also 
explain the Triangulum-A ndromeda structure presented in 
iRocha-Pinto et all ||2003). Yet, given the lack of definitive 
proof of a link between the two structures, we chose not to 
try to fit the radial velocity of the group of Mon ring stars 
behind CMa. Indeed, a small shift of ^ lOkms”^ in the tan¬ 
gential velocities of our simulation substantially changes the 
position and kinematics of the tidal arms. Only with a de¬ 
termination of the radial velocity of the arms away from the 
main body of the dwarf could such a constraint be useful. 


6 CONCLUSION 

We have shown that the Red Giant Branch and Red Clump 
stars at the centre of the CMa dwarf have peculiar kinemat¬ 
ics that are incompatible with a disc-like population: 

(i) The phase space distribution of stars (whether RGB 
or RG) in the CMa fields under the Galactic disc (Figures 3 
and 6) is very different from that of the symmetric region 
on the other side of the disc (Figure 8). 

(ii) More than one third of the stars in the RC sam¬ 
ple follow a distance - radial velocity relation over all the 
sample (5 kpc < Dq < 8 kpc). This relation seems also to 
be present in the RGB sample and prolonged to higher dis¬ 
tances by a group of RGB stars. The important radial ve¬ 
locity shift over only a few klloparsecs for these stars (which 
is not reproduced by the Besangon Galactic model) cannot 
be explained by a group of stars rotating around the Milky 
Way on a circular orbit. 

(iii) The proper motions of the dwarf that we de- 
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rived from the UCAC2.0 catalogue (/iij/rt) = (—3.6 ± 
0.8mas/yr, 1.5 ± 0.4mas/yr) yields at tangential velocities 
when corrected from the Solar motion of {vi,Vb) = (—235 ± 
35kms“^, +15 ± 25kms“^). Combined with a mean radial 
velocity of Vr ~ — 105kms“^ (also corrected from the So¬ 
lar motion), the total velocity of the CMa population is 
Vtot ~ 260kms~^ which is high compared to expecta tions 
for disc stars fe.g. ISoubiran. Bienavme fc Siebertll2003[) . 

(iv) Both the RGB and RC samples show the CMa pop¬ 
ulation is kinematically cold, with intrinsic dispersions of 11 
and 16 kms“^ respectively but this latter value may be over¬ 
estimated due to the difficulty of separating CMa stars from 
contaminating stars in the RC sample. Such low dispersions 
are typical of recent accretion process and can be observed, 
for in stance, in the Sagittarius dwarf galaxy lllhata, et alJ 

liggiD . 

(v) Following the discovery o f a distance spiral arm 
at the edge of the Galactic disc bv lMcClure-GrifEths et mI 
ll2nn4l . we mentioned in Paper III, that the radial velocity 
of the Canis Major RGB population could be compatible 
with this spiral arm. However, the distance - radial veloc¬ 
ity of RG CMa stars reach a much lower value of 72kms~^ 
hardly compatible with the iMcClure-GrifEths et al.l ll2004^ 
values. Moreover, Hi maps do not show any hint of a spi- 
ral arm at the distance of the Canis Major object (see e.g. 
iNakanishi fc Sofuell2n0.'jl . 

All these observations point at the accretion scenario 
for the CMa population, with a dwarf galaxy that is cur¬ 
rently being absorbed by the Milky Way. The low veloc¬ 
ity dispersion suggests that, even though it is undergoing 
dramatic tidal stripping by our Galaxy, the dwarf is still 
bound. But can the CMa dwarf be also responsible for the 
Monoceros ring and the Triangulum-Andromeda structure? 
If simulations can reproduce the observed kinematic signa¬ 
ture of CMa, constraints on the tidal arms of the dwarf are 
at the moment not strong enough to definitely conclude on 
the CMa dwarf being the progenitor of the Mon ring. 

A definite conclusion may only come with future obser¬ 
vations which should concentrate on mapping the regions 
around the centre of dwarf (220° ^ I ^ 260° and —20° < 
b < 0°). Moreover, in addition to constraining the orbit 
of the accretion, determining the intrinsic dispersion of the 
CMa population and the evolution of this dispersion around 
its core should allow a much more precise estimate of the 
parameters of the initial dwarf, as well as an estimate o f its 
mass-loss rate ijohnston. Sigurdsson fc Hernau'istlll99flll . 
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